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Abstract. In this paper, the effect of a centrally located cutout (circular and elliptical) and cracks 
emanating from the cutout on the free fiexural vibration behaviour of functionally graded material 
plates in thermal environment is studied. The discontinuity surface is represented independent of the 
mesh by exploiting the partition of unity method framework. A Heaviside function is used to capture 
the jump in the displacement across the discontinuity surface and asymptotic branch functions are 
used to capture the singularity around the crack tip. An enriched shear flexible 4-noded quadrilateral 
element is used for the spatial discretization. The properties are assumed to vary only in the thickness 
direction. The effective properties of the functionally graded material are estimated using the Mori- 
Tanaka homogenization scheme and the plate kinematics is based on the first order shear deformation 
theory. The influence of the plate geometry, the geometry of the cutout, the crack length, the thermal 
gradient and the boundary conditions on the free fiexural vibration is numerically studied. 

Introduction 

In recent years, a new class of engineered material, the functionally graded materials (FGMs) has 
attracted researchers to investigate its structural behaviour. The FGMs are microscopically inhomo- 
geneous and the mechanical and the thermal properties vary smoothly and continuously from one 
surface to another. FGMs combine the best properties of its constituents. Typically, the FGMs are 
made from a mixture of ceramic and metal. The ceramic constituent provides thermal stability due 
to its low thermal conductivity, whilst the metallic constituent provides structural stability. FGMs 
eliminate the sharp interfaces existing in laminated composites with a gradient interface and are con- 
sidered to be an alternative in many engineering sectors such as the aerospace industry, biomechanics 
industry, nuclear industry, tribology, optoelectronics and other high performance applications where 
the structural member is exposed to high thermal gradient in addition to mechanical load. 

The static and the dynamic characteristics have been studied in detail by many researchers us- 
ing different plate theories, for example, first order shear deformation theory (FSDT) [|29l l38l l34l . 
second and other higher order accurate theory [|24l |6l |2I1 have been used to describe the plate kine- 
matics. Existing approaches in the literature to study plate and shell structures made up of FGMs 
uses finite element method (FEM) based on Lagrange basis functions [|29l |71 |34l, meshfree meth- 
ods [|24l l6l and recently Valizadeh et ah, [|35i used non-uniform rational B-splines based FEM to study 
the static and the dynamic characteristics of FGM plates in thermal environment. Akbari et ah, [|26ll 
studied two-dimensional wave propagation in functionally graded solids using the meshless local 
Petrov-Galerkin method. Huang et al, (H proposed solutions for the free vibration of side-cracked 



FGM thick plates based on Reddy's third-order shear deformation theory using Ritz technique. Ki- 
tipomchai et ah, [fT4l studied nonlinear vibration of edge cracked functionally graded Timoshenko 
beams using Ritz method. Yang et al, [|371 studied the nonlinear dynamic response of a functionally 
graded plate with a through- width crack based on Reddy's third-order shear deformation theory using 
a Galerkin method. Natarajan et ah, [fTSl [T7| and Baiz et ah, O studied the influence of the crack 
length on the free flexural vibrations of FGM plates using the XFEM and smoothed XFEM, respec- 
tively. The above list is no way comprehensive and interested readers are referred to the literature and 
references therein and a recent review paper by Jha and Kant [ 13 1 on FGM plates. 

Plates with cutouts are extensively used in transport vehicle structures. Cutouts are made to lighten 
the structure, for ventilation, to provide accessibility to other parts of the structures and for altering 
the resonant frequency. Therefore, the natural frequencies of plates with cutouts are of considerable 
interest to designers of such structures. Most of the earlier investigations on plates with cutouts have 
been confined to isotropic plates [|22l [Tl[T0l and laminated composites [|28ll32ll . Recently, Janghorban 
and Zare [12J studied the influence of cutout on the fundamental frequency of FGM plates in thermal 
environment using FEM. Their study was restricted to a limited number of configurations, because the 
mesh has to conform to the geometry. In this study, we present a framework that provides flexibility 
to handle internal discontinuities. 

In this paper, we study the influence of a centrally located cutout and cracks emanating from 
the cutouts on the natural frequencies of FGM plates in thermal environment. Circular and ellipti- 
cal cutouts are considered for the study. A structured quadrilateral mesh is used and the cutouts are 
modelled independent of the mesh within the extended finite element (XFEM) framework. A system- 
atic parametric study is carried out to bring the effect of the gradient index {k), the thermal gradient 
(AT = Tc — Tm), the geometry of the cutout on the free flexural vibration behaviour of FGM plates 
with different boundary conditions. 



Theoretical Background 

Reissner-Mindlin plate theory The Reissner-Mindlin plate theory, also known as the first order 
shear deformation plate theory (FSDT) is an extension of the classical plate theory (or, the Kirchhoff 
- Love plate theory). The FSDT takes into account the shear deformations through the thickness and 
are intended for thick plates in which the normal to the medium surface remains straight but not 
necessarily perpendicular to the medium surface. Using the Mindlin formulation, the displacements 
M, f , w at a point (x, z) in the plate (see Figure ([T])) from the medium surface are expressed as 
functions of the mid-plane displacements Mq, fo, Wo and independent rotations O^, 9y of the normal in 
yz and xz planes, respectively, as 

u{x, y, z, t) = Uoix, y, t) + ze^{x, y, t) 
f (x, y, z, t) = Voix, y, t) + zOyi^x, y, t) 

w{x,y,z,t) = Woix,y,t) (1) 
where t is the time. The strains in terms of mid-plane deformation can be written as: 




The midplane strains Ep, the bending strains £b and the shear strains £s in Equation ([2]) are written as: 



a 

Fig. 1 : Coordinate system of a rectangular plate. 
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where the subscript 'comma' represents the partial derivative with respect to the spatial coordinate 
succeeding it. The membrane stress resultants N and the bending stress resultants M can be related 
to the membrane strains, Sp and the bending strains through the following constitutive relations: 
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where the matrices A = Aij, B = Bij and Db = Dij; (i, j = 1,2,6) are the extensional, the bending- 
extensional coupling and the bending stiffness coefficients and are defined as: 
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{Aij, Bij, Dij} = j Q,j{l, z, z^] dz 



(5) 
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The thermal stress resultants, and the moment resultants M* are: 
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where the thermal coefficient of expansion a{z, T) is given by Equation ( 19) and AT(2;) = T{z) — 



To is the temperature rise from the reference temperature To at which there are no thermal strains. 
Similarly, the transverse shear force Q = {Qxz, Qyz] is related to the transverse shear strains Sg 
through the following equation: 



(7) 



where E = Eij = /_(y2 Qij'^t'^j dz; = 4, 5) is the transverse shear stiffness coefficient, Vi, Vj is 
the transverse shear coefficient for non-uniform shear strain distribution through the plate thickness. 
The stiffness coefficients Q^, are defined as: 
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(8) 



where the modulus of elasticity E{z, T) and Poisson's ratio u are given by Equation (18). The strain 
energy function U is given by: 



where 6 = {u, v,w,9x,9y} is the vector of the degrees of freedom associated to the displacement field 
in a finite element discretization. Following the procedure given in [i27il . the strain energy function U 
given in Equation ([9]) can be rewritten as: 

U{S) = ^<5^K<5 (10) 
where K is the linear stiffness matrix. The kinetic energy of the plate is given by: 

T{S) = \j {p{ul + vl + wl) + i{el + el)] d^i (11) 

where p = J^f^^^ P^^) ^ ^ /^h/2 ^^'^ P(^) '^^^ mass density that varies through 

the thickness of the plate. The external work due to the in-plane stress resultants {N^^^i -^to' 
developed in the plate under the thermal load is: 



V{6) 



(12) 



Upon substituting Equations (10) - (12) in Lagrange's equations of motion, the following governing 
equation is obtained: 



+ (K + Kg)5 = 



(13) 



where M is the consistent mass matrix. After substituting the characteristic of the time function [8] 
d = —co'^S, the following algebraic equation is obtained: 



((K + Kg) - w^M) 6 = 



(14) 



where K and Kg are the stiffness matrix and the geometric stiffness matrix due to thermal loads, 
respectively, u is the natural frequency. The plate is subjected to a temperature field. So, the first step 



in the solution process is to compute the in-plane stress resultants due to the temperature field. These 
will then be used to compute the stiffness matrix of the system and then the frequencies are computed 
for the system. 

Functionally Graded Material A functionally graded material (FGM) rectangular plate (length 
a, width b and thickness h), made by mixing two distinct material phases: a metal and a ceramic is 
considered with coordinates x, y along the in-plane directions and z along the thickness direction (see 
Figure ([T])). The material on the top surface (z = +h/2) of the plate is ceramic and is graded to metal 
at the bottom surface of the plate (z = —h/2) by a power law distribution. The homogenized material 
properties are computed using the Mori-Tanaka Scheme [fT5l [3l. 

Estimation of mechanical and thermal properties Based on the Mori-Tanaka homogenization 
method, the effective bulk modulus K and shear modulus G of the FGM are evaluated as IIT5l l3ll5ll25l 



where 




Fig. 2: Through thickness variation of volume fraction 



Here, Vi {i = c, m) is the volume fraction of the phase material. The subscripts c and m refer to 
the ceramic and the metal phases, respectively. The volume fractions of the ceramic and the metal 
phases are related by \4 + = 1 and Vc is expressed as 

K(^)=(^)', fc>0 (17) 



where k in Equation (17) is the volume fraction exponent, also referred to as the gradient index. 
Figure (|2]) shows the variation of the volume fraction of ceramic constituent in the thickness direction 
z for a FGM plate. The effective Young's modulus E and Poisson's ratio u can be computed from the 
following expressions: 

9KG _ 3K-2G 

^-WTg' ''-2{m + Gy ^^^^ 



The effective mass density p is given by the rule of mixtures as: p = p^Vc + PmVm- The effective heat 
conductivity k^s and the effective coefficient of thermal expansion a^s is given by: 
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(19) 



where = k^- Km and Ctcm = ftc - ttr 



Temperature distribution through the thickness The material properties P that are temperature 
dependent can be written as [30|: 



P = Po{P-iT-^ + 1 + PiT + PsT' + PsT^), 



(20) 



where Po, P-i, Pi, P2, P3 are the coefficients of temperature T and are unique to each constituent 
material phase. The temperature variation is assumed to occur only in the thickness direction and 
the temperature field is considered to be constant in the xy-plane. In such a case, the temperature 
distribution along the thickness can be obtained by solving a steady state heat transfer problem: 
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The solution of Equation (21 1 is obtained by means of a polynomial series [|36ll as: 

T{z) =Tm + {Tc - Tm)r]{z, h) 



(21) 
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Element Formulation 



Shear flexible Q4 element The plate element employed here is a C° continuous shear flexible field 
consistent element with five degrees of freedom (mq, Vo,Wo,0x,6y) at four nodes in a 4-noded quadri- 
lateral (QUAD-4) element. If the interpolation functions for QUAD-4 are used directly to interpolate 
the five variables (uo, Vo, Wo, 9x, Oy) in deriving the shear strains and the membrane strains, the ele- 
ment will lock and show oscillations in the shear and the membrane stresses. The field consistency 
requires that the transverse shear strains and membrane strains must be interpolated in a consistent 



manner. Thus, the 6^ and 6y terms in the expressions for shear strain Ss have to be consistent with the 
derivative of the field functions, Wo^x and Wo,y. This is achieved by using field redistributed substitute 
shape functions to interpolate those specific terms, which must be consistent as described in [[33l [8l. 
This element is free from locking and has good convergence properties. For complete description of 
the element, interested readers are referred to the literature [j33l[8^. where the element behaviour is dis- 
cussed in great detail. Since the element is based on the field consistency approach, exact integration 
is applied for calculating various strain energy terms. 

Enriched Q4 element In this study, an XFEM framework is employed to represent the internal 
discontinuity. XFEM is classified as one of the partition of unity methods. Using this property, any 
function Lp can be reproduced by a product of the partition of unity shape functions with (/?. Let 
u'^ C U, the XFEM approximation can be decomposed into the standard part u^^^ and into an enriched 
part u^„, as: 

u'^(x) = ut,(x) + u,V(x) 

= E ^i(x)^(x)aj (25) 

where A/''^^™ is the set of all the nodes in the EE mesh, M'^™ is the set of nodes that are enriched with 
the enrichment function g. The enrichment function g carries with it the nature of the solution or the 
information about the underlying physics of the problem, for example, g = H, is used to capture 
strong discontinuities, where H is the Heaviside function or g could be a set of functions that are 
chosen to represent the near tip asymptotic fields. The necessary steps involved in the implementation 
of the XFEM are: 

1. Representation of the interface In this study, a level set approach is followed to model the 
cutouts. The geometric interface (for example, the boundary of the cutout) is represented by 
the zero level curve = 0(x, t) = 0. The interface is located from the value of the level set 
information stored at the nodes. The standard EE shape functions can be used to interpolate 
at any point x in the domain as: 

0(x) = EiV,(x)0, (26) 

where the summation is over all the nodes in the connectivity of the elements that contact x 
and (pj are the nodal values of the level set function. For circular cutout, the level set function 
is given by: 

0/ = ||x/ - Xcll - rc (27) 

where Xc and are the center and the radius of the cutout. For an elliptical cutout oriented at 
an angle 9, measured from the x— axis the level set function is given by: 

(pi = \/ai{xi - Xc)2 _ a2{xi - Xc){yi - Vc) + a^ivi - Vc)^ - 1 (28) 

where 



where d and e are the major and minor axes of the ellipse and {x^, Hc) is the center of the ellipse. 
Figure (|3]) shows the level set function for a circular and an elliptical cutout. 




(b) Elliptical cutout 
Fig. 3: Level set information for different shapes. 



2. Choice of enrichment function In this study, the enriched approximation takes the following 
form: 
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where A/"^*^™ is the set of all nodes in the finite element mesh, M^.M^ and AT'^"* are the set of 
nodes enriched with a Heaviside function, the near-tip asymptotic fields and the absolute value 
function, aj and are the nodal degrees of freedom corresponding to the Heaviside function 
H and the near-tip functions, {B.a]i<a<i, cl is the nodal degrees of freedom that corresponds 
to the enrichment function ^'(x) used to represent the jump in the displacement to represent the 
cutouts. In this study, a Heaviside function is used to capture the physics of the problem (viz., 
jump across the discontinuity surface) and a set of asymptotic functions that represent the near 
tip fields are used, given by: 

{ 6 6 6 6^ 

{E:„}i<„<4(r, 61) = <^ sin -, cos -, sin 6 am.-, sin 6' cos - > (32) 

3. Integration The standard GauB quadrature cannot be applied in elements enriched by discon- 
tinuous terms, because GauB quadrature implicitly assumes a polynomial approximation. In the 
present study, a triangular quadrature with sub-division aligned to the discontinuous surface is 
employed. For the elements that are not enriched, a standard 2x2 Gaussian quadrature rule 
is used. The other techniques that can be employed are Schwarz Christoffel Mapping [fT9ll20l , 
Generalized quadrature [[T6l and Smoothed extended FEM [4J. 



Remark. Appropriates terms from the displacement approximation, given by Equation (30) are cho 



sen for the study. For example, the first and the last terms are retained for the problems with cutouts 



only, whilst, all the terms in Equation (30) are used for the problem with cracks and cutouts. 



Remark. The degrees of freedom of the unused nodes in the EE mesh are taken care of during the 
solution process. 



Numerical Results 

In this section, we study the influence of a centrally located cutout on the fundamental frequencies 
of FGM plates. We consider both square and rectangular plates with simply supported and clamped 
boundary conditions. Two different cutout shapes, viz., circular and elliptical cutouts (see Figure (|4])) 
are considered in this study. Table [T] presents the integration rules used for the current study. In all 
cases, we present the non dimensionalized free flexural frequencies as, unless specified otherwise: 



(a) 



(b) 




X 



Fig. 4: Plate with a centrally located circular and an elliptical cutout, r is the radius of the circular 
cutout, 2d and 2e are the major and minor axes defining the ellipse. 



Table 1: Integration rules for enriched and non-enriched elements. 



Element Type 



Non-enriched element 

Tip element 

Tip blending element 

Split element 

Split blending element 



GauB points 



13 per triangle 
16 

3 per triangle 
4 



Split- Tip blending element 4 per triangle 



Table 2: Temperature dependent coefficient for material Si3N4/SUS304, Ref [|30ll34l. 



Material Property Po P-i Pi P2 P3 





E(Pa) 


348.43e9 


0.0 


-3.070e-* 


2.160e-^ 


-8.946e-i^ 


Si3N4 














a (1/K) 


5.8723e-6 


0.0 


9.095e-4 


0.0 


0.0 




E(Pa) 


201.04e9 


0.0 


3.079e-^ 


-6.534e-^ 


0.0 


SUS304 














a (1/K) 


12.330e-*^ 


0.0 


8.086e-^ 


0.0 


0.0 



where Dc = Ech^ / (12(1 — u"^)) and pc is the mass density. In order to be consistent with the existing 
literature, properties of the ceramic are used for normalization. The effect of the plate slendemess 
ratio a/h, the plate aspect ratio h/a, the cutout radius r/a, the cutout geometry d/e and the boundary 
condition on the natural frequencies are numerically studied. The FGM plate considered here consists 
of silicon nitride (Si3N4) and stainless steel (SUS304). The material is considered to be temperature 
dependent and the temperature coefficients corresponding to Si3N4/SUS304 are listed in Table |2] [34, 
1301 . The mass density (p) and the thermal conductivity {K) are: pc = 2370 kg/m^, Kc = 9.19 W/mK 
for Si3N4 and pm = 8166 kg/m^, = 12.04 W/mK for SUS304. Poisson's ratio u is assumed to be 
constant and taken as 0.28 for the current study [|34l l23ll . Here, the modified shear correction factor 
obtained based on energy equivalence principle as outlined in [1311 is used. The boundary conditions 
for simply supported and clamped cases are : 
Simply supported boundary condition: 

Uo = Wo = Oy = on X = 0, a; Vo = Wo = 9x = on y = 0,b (34) 
Clamped boundary condition: 

Uo = Wo = 9y = Vo = 0^ = on X = 0,a Sz y = 0,b (35) 



Validation Before proceeding with a detailed study on the effect of different parameters on the 
natural frequency, the formulation developed herein is validated against available results pertaining 
to the linear frequencies of a FGM plate in thermal environment [11 J and an isotropic plate with 
a centrally located circular and an elliptical cutout [|10]. The computed linear frequencies: (a) for a 
square simply supported FGM plate in thermal environment with a/h = Sis given in Table |3] and (b) 
for square plate with a circular cutout for various boundary conditions is given in Table |4j It can be 
seen that the numerical results from the present formulation are found to be in good agreement with 
the existing solutions. Based on a progressive refinement a 40 x 40 quadrilateral mesh is found to be 
adequate to model the full plate with a cutout. 



Effect of thermal environment Next, the linear free flexural vibration behaviour of FGM is numer- 
ically studied with and without thermal environment. For the uniform temperature case, the material 
properties are evaluated at Tc = Tm = 300K. The temperature field is assumed to vary only in 
the thickness direction and determined by Equation ([22]). The temperature for the ceramic surface 



Table 3: Comparison of non-dimensional linear frequency ^ - 
ported FGM plate {a/b = 1, a/h = 8) in thermal environment. 
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Temperature 


gradient index 


Mode 1 


Mode 2 




T T 


k 


Ref. [11] 


Present 


Ref. [llj 


Present 






0.0 


12.397 


12.315 


29.083 


29.031 


T 


= 400K, Tm = 300K 


0.5 


8.615 


8.484 


20.215 


19.986 


1.0 


7.474 


7.443 


17.607 


17.515 






2.0 


6.693 


6.679 


15.762 


15.709 






0.0 


11.984 


11.894 


28.504 


28.436 


T 


= 600K, Tm = 300K 


0.5 
1.0 


8.269 
7.171 


8.147 
7.126 


19.784 
17.213 


19.535 
17.098 






2.0 


6.398 


6.370 


15.384 


15.309 



Table 4: Natural frequency parameter 
cutout. 



n 



Cutout 



Circular (r /a = 0.1) 



1/4 



for an isotropic plate with a central 



Boundary 
Condition 



Mesh 



Frequency 



mode 1 mode 2 



cccc 



Elliptical, (2d/a =3/8; 2d/2e = 2) SSSS 



20x20 
30x30 
40x40 

Ref. niOl 
20 x20 
30x30 
40x40 

Ref. ttlOl 



6.1848 
6.1762 
6.1725 
6.2110 
4.4828 
4.4775 
4.4758 
4.4820 



8.7215 
8.6622 
8.6443 
8.7310 
6.9237 
6.8849 
6.8705 
6.9120 



is varied (Tc = 400 K, 600 K, 900 K), whilst maintaining a constant value on the metallic surface 
{Tm = 300K) to subject the plate to a thermal gradient. Figure ^ shows the transverse displacement 
along y = 6/2 for a simply supported square FGM plate under different thermal loading with center 
circular cutout r/a = 0.2. It can be seen that the transverse displacement increases with increasing 
thermal gradient. The transverse displacement for a FGM plate with gradient index A; = 2 is greater 
than a plate with k = (note that = is a pure ceramic plate) as expected. Figure (|6]) shows the 
influence of cutout radius on the maximum transverse deflection for a simply supported square FGM 
plate in thermal environment, AT = lOOK. It can be seen that increasing the gradient index, increases 
the transverse deflection. With increase in the cutout radius, the maximum transverse displacement 
first increases due to stiffness reduction and with further increase the transverse displacement de- 
creases. This is because for a larger cutout radius, there is less material to deform. The increase in 
the transverse displacement with increasing gradient index can be attributed to the increase in the 
metallic volume fraction, which increases the coefficient of thermal expansion and thus increasing 
the transverse displacement. The increase in the displacement with cutout radius is due to stiffness 
degradation. 
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Fig. 5: Transverse displacement along the centerline {w{x, b/2, 0)) of the plate for various gradient 
indices (k = 0, 2) and for different thermal gradient, viz., AT = OK, lOOK, 300K (T^ = 300K, 400K 
and 600K and T^ = 300K) with a circular cutout (r/a = 0.2) at the center. 

The geometric stiffness matrix is computed from the in-plane stress resultants due to the applied 
thermal gradient. The geometric stiffness matrix is then added to the stiffness matrix and then the 
eigenvalue problem is solved. Table |5] shows the influence of the gradient index, the thermal gradient 
and the plate aspect ratio on the fundamental frequency of a square simply supported FGM plate with 
a circular cutout r/a = 0.2. It can be seen that the combined effect of increasing the gradient index 
and the thermal gradient is to lower the fundamental frequency, whilst the frequency increases with 
decreasing plate thickness. The effect of boundary conditions, the plate thickness and the gradient 
index k on the fundamental frequency for a square FGM plate in thermal environment (AT = lOOK) 
is shown in Table [6j The fundamental frequency decreases with increasing gradient index due to 




Fig. 6: Maximum transverse displacement {wmax = w/h) as a function of the cutout radius for a 
square simply supported FGM plate for different gradient index {k = 0, 2) subjected to a thermal 
gradient AT = WOK. 

increase in the metallic volume fraction. The frequency initially increases with plate thickness but 
upon further increase, the fundamental frequency decreases. Figure (|7]) shows the influence of the 
plate aspect ratio on the fundamental frequency for a FGM plate in thermal environment AT = lOOK 
for different boundary conditions with gradient index k = 2 and a/h = 10. It can be seen that the 
frequency increases with increasing a/b ratio and clamped plate has higher frequency than a simply 
supported plate. The increase in the stiffness is the cause for increase in frequency when the boundary 
condition is changed from simply supported to clamped condition for a fixed aspect ratio and plate 
thickness. 



Effect of cutout geometry Figure ([8]) shows the influence of the cutout size on the frequency for 
a plate in thermal environment (AT = lOOK). The frequency increases with increasing cutout size, 
whilst decreases with increasing gradient index. The effect of geometry of the cutout d/eis shown in 
Figure (|9]) for a square FGM plate in two different thermal environment, viz., AT = OK, lOOi^. It 
can be seen that the fundamental frequency increases with increasing the d/e ratio, whilst decreases 
with increasing thermal gradient. 

Table|7] shows the influence of the orientation 9 of an elliptical cutout on the fundamental frequen- 
cies for a square simply supported FGM plate in thermal environment with a/h = 10 and gradient 
index k = 2. The geometry of the plate is shown in Figure (10). It can be seen that the orientation 9 
of the cutout has strong influence on the fundamental frequencies. It can be seen that increasing the 
orientation 9 from 0° to 90°, the mode 1 and the mode 3 first decreases until 9 = 45° and with further 
increase in 9, the mode 1 and mode 3 frequency increases. On the other hand, the mode 2 frequency 
first increases until 9 = 45° and then decreases. The frequencies are symmetric with respect to 6* = 
45°, indicating that the FGM plate is globally homogeneous. With increasing temperature, all the 
fundamental frequencies decreases as observed earlier. 



Table 5: Non-dimensionalized mode 1 frequency for square plates with circular cutouts r/a = 
with a/h — 5,\Q 'm thermal enviromnent. 



a/h 


k 


Tc = 300K 


= 400K 


% = 600K 


% = 900K 







17.6855 


17.4690 


17.0266 


16.3111 




1 


10.6681 


10.5174 


10.1932 


9.6350 


5 


2 


9.6040 


9.4618 


9.1469 


8.5882 




5 


8.7113 


8.5738 


8.2544 


7.6601 




10 


8.2850 


8.1484 


7.8191 


7.1840 







19.1844 


18.5992 


17.2928 


14.8701 




1 


11.5736 


11.1317 


10.1161 


8.1589 


10 


2 


10.4135 


9.9844 


8.9842 


7.0277 




5 


9.4461 


9.0145 


7.9853 


5.9172 




10 


8.9858 


8.5452 


7.4747 


5.2682 



Table 6: Non-dimensionalized mode 1 linear frequency for square plates with circular cutouts r/( 
0.2, temperat ure gradient AT = 100 K (7; = 400K, = 300K). 

Boundary a/h gradient index. A; 



Condition 







1 


2 


5 


10 




5 


17.4690 


10.5174 


9.4618 


8.5738 


8.1484 


SSSS 


10 


18.5992 


11.1317 


9.9844 


9.0145 


8.5452 




20 


17.5380 


10.1776 


9.0071 


7.9664 


7.4294 




25 


16.3587 


9.1762 


7.9919 


6.8810 


6.2689 




5 


31.4944 


18.9259 


16.9575 


15.3461 


14.6221 


cccc 


10 


38.7777 


23.2411 


20.8571 


18.9072 


18.0010 




20 


41.0541 


24.4016 


21.8525 


19.7479 


18.7389 




25 


40.7846 


24.0796 


21.5110 


19.3700 


18.3275 




Fig. 7: Effect of the plate aspect ratio on the linear frequency for a FGM plate with a/h =10 and 
with a thermal gradient AT = lOOK = 400K, T„ = 300K) for a centrally located circular cutout 
with radius r/a = 0.2 for various boundary conditions (simply supported and clamped) with gradient 
index k = 2. Note that in this study, a is kept constant. 




Fig. 8: Effect of the cutout size on the linear frequency (Q) for a square simply supported FGM plate 
with a/h =10 in a thermal environment AT = lOOK (Tc = 400K, T^ = 300K) for a centrally located 
cutout for various gradient indices k. 
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Fig. 9: Influence of the geometry of the cutout on the linear frequency {fl) for a square FGM plate 
with a/h =10 with gradient index /c = 2 for various thermal gradients. 




Fig. 10: Plate with a centrally located elliptical cutout oriented at an angle 6, measured in the coun- 
terwise direction from the x— axis. 



Table 7: Non-dimensionalized natural frequency for a square simply supported FGM plate (a/6 = 1, 
a/h = 10, /c = 2) with an elliptical cutout. 

Major axis = 300K = 400K 



orientation mode 1 mode 2 mode 3 mode 1 mode 2 mode 3 



0° 


9 


6795 


16 


3035 


23 


0825 


9 


2771 


16 


0082 


22 


7147 


10° 


9 


6627 


16 


3067 


23 


0677 


9 


2591 


16 


0108 


22 


6994 


20° 


9 


6231 


16 


3311 


23 


0329 


9 


2168 


16 


0338 


22 


6632 


30° 


9 


5740 


16 


3353 


22 


9892 


9 


1643 


16 


0364 


22 


6175 


40° 


9 


5435 


16 


3548 


22 


9596 


9 


1317 


16 


0548 


22 


5868 


50° 


9 


5435 


16 


3548 


22 


9596 


9 


1317 


16 


0548 


22 


5868 


60° 


9 


5740 


16 


3353 


22 


9892 


9 


1643 


16 


0364 


22 


6175 


70° 


9 


6231 


16 


3311 


23 


0329 


9 


2168 


16 


0338 


22 


6632 


80° 


9 


6627 


16 


3067 


23 


0677 


9 


2591 


16 


0108 


22 


6994 


90° 


9 


6795 


16 


3035 


23 


0825 


9 


2771 


16 


0082 


22 


7147 



Cracks emanating from the cutout Next, we study the influence of cracks emanating from an 
elliptical cutout. Two cases are considered as show in Figure (fTTl), viz., (a) Case A: elliptical cutout 
with two cracks one on either side with crack length i = 0.25a (b) Case B: single elliptical cutout 
with length of the major axis equals the sum of the lengths of the two cracks and an elliptical cutout 
in Case A. Table [8] shows the influence of the gradient index and the thermal gradient on the natural 
frequencies of a square simply supported FGM plates. It can be seen that increasing the temperature 
of the ceramic phase and the gradient index, the fundamental frequency decreases. The mode 1 and 
mode 3 frequency for both the cases show very similar behaviour, whilst the mode 2 frequency for 
Case A is greater than the mode 2 frequency for Case B. This can be attributed to the fact that Case A 
can considerably more material to resist the deformation and that in the case of mode 2, the crack is 
parallel to the deformed shape. The mode 2 deformed shape can be seen as the crack opening mode. 
On the other hand, the crack and the cutout is perpendicular to the deformed shape in the case of 
mode 3, which resists the crack opening. 



Next, the influence of multiple cracks emanating from a circular cutout is studied. Figure (12) 
shows the influence of the crack orientation on the first fundamental frequency. In this case, a simply 
supported square FGM plate with gradient index. A; = 2, a centrally located circular cutout r/a = 0.2 
and crack with length £/a = 0.1 is considered. It can be seen from Figure (12) that the fundamental 
frequency is symmetric with respect to crack orientation 9 = 45° and has a least value at 6* = 45°. 
Figure ( 13 ) shows the influence of number of cracks emanating from a circular cutout on the funda- 
mental frequency for a simply supported square FGM plate with a/h = 10, gradient index k = 2, 
cutout ratio r/a = 0.2 and exposed to different thermal gradients. It can be seen that with increasing 
number of cracks and thermal gradient, the fundamental frequency decreases. This can be attributed 
to the increase in the local flexibility due to the presence of cracks and thus decreasing the frequency. 



0.1a 



0.1a 



(a) (b) 

Fig. 11: Simply supported FGM plates a//i = 10, a) with an elliptic cutout and two cracks b) with an 
elliptic crack. 

Table 8: Non-dimensionalized natural frequency for square FGM plates with: (a) elliptic cutout with 
two cracks and (b) elliptic crack 

Case A: ellipse with two cracks Case B: elliptic crack 

AT k mode 1 mode 2 mode 3 mode 1 mode 2 mode 3 

16.2790 34.4610 41.5259 16.1708 20.4605 42.7517 

2 8.8710 18.6648 22.5162 8.8104 11.1051 23.1853 

5 8.0473 16.8857 20.4443 7.9928 10.0659 21.0584 

15.7865 34.0478 40.9831 15.6527 20.0172 42.1728 

100 2 8.5066 18.3994 22.1638 8.4238 10.7881 22.8061 

5 7.6827 16.6294 20.0996 7.6057 9.7508 20.6869 



Influence of the crack length As a last example, the influence of the crack length on the first three 
fundamental frequency is studied. The geometry of the cutout, the location of the crack and the plate 
geometry is shown in Figure ([14]). The plate is exposed to a thermal gradient, AT = 300K (T^ = 
600K and T^ = 300K). It can be seen from Figure ( 14) that increasing the crack length decreases the 
fundamental frequencies. The decrease in the mode 1 frequency is linear whilst the mode 3 frequency 
decreases in a bi-linear mode. For the rectangular plate considered here, the effect of increasing crack 
length has minimal impact on the mode 2 frequency. This is because for the rectangular plate (see 
Figure ([14]), the crack is perpendicular to the deformed shape. 



Summary 

The influence of a centrally located cutout and cracks emanating from the periphery of the cutout on 
the fundamental frequency of FGM plates is numerically studied. The influence of various parameters, 
viz., the plate thickness h, the plate aspect ratio a/b, the cutout geometry d/e, the cutout radius 
r/a, the crack length i/a and the gradient index k is studied for FGM plates in different thermal 
environment. The formulation is based on the first order shear deformation theory for plates and an 
enriched four-noded field consistent shear flexible element is used. The material is assumed to be 



10.4 




10 20 30 40 50 60 70 80 90 

e (degrees) 

Fig. 12: The effect of the crack orientation 9 on the mode 1 natural frequency (Vl) for a square simply 
supported FGM plate (a/h =10, k =2) with a circular cutout r/a = 0.2 and AT =0 (Tc = Tm = 
300K). The length of each crack is £/a = 0.1. 
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7.5 1 ' ' ' • ' 1 
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Fig. 13: The effect of the number of cracks on the mode 1 natural frequency (fi) for a square sim- 
ply supported FGM plate (a/h = 10, k =2) with a circular cutout r/a = 0.2 in different thermal 
environments. The length of each crack is £/a = 0.1. 

temperature dependent and graded in the thickness direction. From a systematic parametric study, the 
following can be concluded: 

• Increasing the cutout radius r/a and the cutout geometry d/e increases the frequency. This can 
be attributed to stiffness degradation. 
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Fig. 14: The effect of the crack length £/a on the first three fundamental frequencies (Q) for a rect- 
angular simply supported FGM plate (b/a — 2, a/h — 10, k — 2) with a circular cutout r/a — 0.2 in 
a thermal environment with AT =300K. 

• Increasing gradient index k decreases the natural frequency. This is due to the increase in the 
metallic volume fraction. 

• Increasing the slendemess ratio a/h and the aspect ratio a/b increases the frequency. 

• Increasing the crack length i/a and increasing the number of cracks, decreases the natural 
frequency. 

• Increasing the thermal gradient decreases the frequency. 

It can be observed that the combined effect of increasing the cutout radius, the cutout geometry 
and decreasing gradient index is to increase in the frequency. In all the cases, the increase is due to 
stiffness degradation. In case of the gradient index, decreasing metallic volume fraction increases the 
stiffness and thus increasing the frequency. 
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